Quantum Trajectories for Realistic Photodetection II: Application and Analysis 
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In the preceding paper [Warszawski and Wiseman] we presented a general formalism for deter- 
mining the state of a quantum system conditional on the output of a realistic detector, including 
effects such as a finite bandwidth and electronic noise. We applied this theory to two sorts of 
photodetectors: avalanche photodiodes and photoreceivers. In this paper we present simulations of 
these realistic quantum trajectories for a cavity QED scenario in order to ascertain how the con- 
ditioned state varies from that obtained with perfect detection. Large differences are found, and 
this is manifest in the average of the conditional purity. Simulation also allows us to comprehen- 
sively investigate how the quality of the the photoreceiver depends upon its physical parameters. In 
particular, we present evidence that in the limit of small electronic noise, the photoreceiver quality 
can be characterized by an effective bandwidth, which depends upon the level of electronic noise 
and the filter bandwidth. We establish this result as an appropriate limit for a simpler, analytically 
solvable, system. We expect this to be a general result in other applications of our theory. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc, 42.50.Ar 



I. INTRODUCTION 

In the preceding paper |jj we gave a full description 
of a method to model the evolution of open quantum 
systems conditional upon detection results from realis- 
tic detectors. This is a generalization of standard quan- 
tum trajectory theory, to take into account correlations 
between the system and classical detector states which 
cannot be observed in practice. We also showed how it 
can be applied in quantum optics, deriving realistic quan- 
tum trajectories for conditioning upon photon counting 
using an avalanche photodiode (APD), and homodync 
detection using a photoreceiver (PR). The greatest sig- 
nificance of our work is in the field of quantum control, 
where conditional states are the optimal basis for control 
loops. 

In this paper we find and study solutions to the re- 
alistic quantum trajectory equations derived in the pre- 
ceding paper. We consider the evolution of a two-level 
cavity QED system (which for convenience we refer to 
as a two- level atom), conditioned on four different types 
of detection (using the two detectors mentioned above). 
Some of these solutions have been studied in a previous 
paper by us and Mabuchi Q , but not with the thorough- 
ness we apply here. These solutions can only be found 
numerically, because of the nonlinearity of the system. 
We also consider another system, the degenerate para- 
metric oscillator below threshold (DPOBT), which can 
be treated analytically for realistic homodyne detection. 

Our study in this paper achieves five important aims. 
First, it shows how the theory developed in the preceding 
paper is implementable in practice. Second, it establishes 
the degree of impact of detection imperfections on the pu- 
rity of the conditional system state. Third, it illustrates 
how realistic quantum trajectories typically differ from 
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standard (idealized) quantum trajectories, which helps 
to build intuition about them. Fourth, it emphasizes 
the importance of the experimenter's choice of detection 
scheme in determining how well one can follow the condi- 
tional evolution of a system. Fifth, it allows us to quan- 
titatively investigate emergent properties of realistic tra- 
jectories, such as the "effective bandwidth" discussed in 
the preceding paper. 

This paper is organized as follows. We begin in Sec. II 
by introducing the two-level atom (TLA) that will be the 
subject of study for most of the remainder of the paper. 
We also consider standard (idealized) quantum trajecto- 
ries for this system under the four different measurement 
schemes we consider. In Sec. Ill we discuss the idea of 
levels of conditioning. These arise from different levels 
of knowledge a hypothetical observer has about the dy- 
namics of the detector. They help one to understand how 
realistic trajectories differ from idealized trajectories. In 
Sec. IV we analyze the stochastic dynamics of the system 
under two realistic detection schemes (direct and adap- 
tive) using the APD. In Sec. V we do likewise for the two 
homodyne detection schemes using the PR. In Sec. VI 
we return to the question of effective bandwidth for the 
PR, discussed in Sec. IV B of the preceding paper. We 
verify the formula suggested in the preceding paper, nu- 
merically for the TLA and analytically for the DPOBT. 
Sec. VII contains a discussion of the numerical techniques 
used in our simulations, and Sec. VIII concludes. 



II. THE TWO-LEVEL SYSTEM 
A. The system 

A driven (0) and damped (7) two-level atom (TLA) 
is the smallest (in Hilbert space) quantum optical sys- 
tem. However, it is a nonlinear system and has surpris- 
ingly rich dynamics. These are especially evident when 
one considers conditioning on measurements of its fluo- 
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rescence ||, |[ |?| |(|. For this reason, we use it as a test 
system for investigating realistic photodetection. 

In reality it would be almost impossible to collect and 
detect all, or even most, of the fluorescence of an atom in 
free space. It is therefore better to envisage our two-level 
quantum system as a cavity- QED system consisting of a 
single two-level atom strongly coupled (<?) to an optical 
cavity which is even more strongly damped (k). In the 
limit k 3> 3, f2 3> 7, the light emitted by the atom into 
the cavity damps through the cavity end-mirror to form 
an output with the same temporal properties as the free- 
space fluorescence. If g 2 / n ^S> 7, the free-space damping 
may be ignored and the effective damping rate of the 
two- level system is T = g 2 / k. We have this in mind when 
choosing parameters for our simulations. The advantage 
of this scheme is of course that the cavity output beam 
is readily measurable, by photon counting or homodyne 
detection. 

Let us denote the ground and excited states for the 
TLA by I g) and |e). Defining the Pauli matrices in the 
usual manner, the state matrix for the TLA can be writ- 
ten as 



p = \ {I + xa x + yijy + za z ) , 



(2.1) 



where (x, y, z) is the Bloch vector which is confined to a 
unit-sphere. The purity of the TLA can be defined as 



purity = Tr[p 2 ] = \ (l 4- x 2 + y 2 + z 2 ) 



(2.2) 



It has an upper limit of 1 (a pure state) and lower limit 
of I (a completely mixed state). 

The unconditioned ME for the driven and damped 
TLA in the interaction picture is 

p=-i-[a x ,p]+V[c]p = Cp. (2.3) 

Here c = \/Ta, where a is the TLA lowering operator. 
The steady state solution of this ME is 
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p = 



&T<t v +r 2 (1 - cr 2 )/2 



which has a purity of 



2n 2 + r 2 



n 2 



2n 2 + r 2 



(2.4) 



(2.5) 



which goes to 1 for small and to i for large tt. 
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FIG. 1: A typical trajectory for direct detection for £1 = 10F. 
The two quadratures of the TLA are shown (x in plot (A) 
and y in plot (B)) as is the inversion, z, of the TLA (in plot 
(C). Time is measured in units of This is true for the 

remainder of this paper, unless otherwise stated. 



Here E[dN(t)} = Tr[( C t + p*)(c + p)p]dt 1 and fi is the 
(necessarily small) local oscillator (LO) amplitude. See 
the preceding paper for definitions of the superoperator 
symbols Ti and Q . A typical quantum trajectory for di- 
rect detection (for which p — 0) is shown in Fig. |l|, where 
plots of x, y, z specify the system state. 

The x quadrature is zero since the Hamiltonian is pro- 
portional to this quadrature. For this reason we can en- 
visage the TLA moving in the y-z plane of the Bloch 
sphere. Since 3> T for this trajectory, the y value of the 
state oscillates between —1 and 1 between jumps. The 
coherent driving also causes the inversion of the TLA, 
z, to oscillate between —1 (ground state) and 1 (excited 
state). Jumps (photon detections) take z to — 1 and y 
to as expected. The conditional purity is of course 
unity at all times. The average phot on f lux entering the 
detector can be calculated from Eq. (2.4) to be 
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2n 2 + l 2 



(2.7) 



This saturates at F/2 for 17 ^ T, as is consistent with 
the number of jumps seen in Fig. H. 



B. Perfect Direct Detection 

Let us now consider continuous, perfect measurement 
of the TLA. Firstly consider counting photons in the 
emitted field. In the preceding paper we stated the rele- 
vant stochastic master equation (SME) in Sec. II C. With 
efficiency r\ = 1 it is 

dp = - dtH[iH + \c x c + p* c + \\p\ 2 ]p 

+ dN(t)g[(c + p)]p. (2.6) 



C. Perfect Adaptive Detection 

A quite different sort of conditional dynamics are re- 
vealed using the adaptive detection scheme proposed in 
H . In this scheme the output field from the TLA is com- 
bined with a weak LO at a low reflectivity beam splitter, 
before being subjected to direct detection (see Fig. ||). 
The local oscillator amplitude p in Eq. ( |2.6| ) is chosen 
to be iVr/2. Two values are present here because the 
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Photodetector 




FIG. 2: Adaptive detection, with the low intensity LO ampli- 
tude controlled by the output of the photodetector. See text 
for further details. 
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FIG. 3: A typical trajectory for adaptive detection with S7 
10r. Two state jumping is clearly visible. 



LO amplitude is flipped following each detection. This 
is what makes the scheme adaptive: the measurement 
being made at any particular time depends on the past 
results. 

With this fi = ±\/r/2 and the Q,a x /2 Hamiltonian, 
the behaviour of the TLA is very simple. After transients 
have passed, the state just jumps between two fixed states 
that are close to x eigenstates (for SI ^> T). These two 
states are actually the stable eigenstates of the two (±) 
no-detection measurement operators (see, for example, 
Ref. [^| ) When a detection on the combined field occurs, 
the TLA jumps into the stable eigenstate of the other 
no jump operator. A typical trajectory for adaptive de- 
tection is shown in Fig. || for SI = 10r. The two-state 
jumping is clearly evident. Note that y and z take on the 
same values in both of the two stable states. 

The actual photon flux entering the photodetector for 
perfect adaptive detection can be calculated using the 
eigenstates of the no-jump operator, which are given in 



FIG. 4: A typical trajectory for homodyne x detection with 

q. = tor. 



Ref. ||. In either eigenstate it is 
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Again, this is consistent with the number of jumps in 
Fig. |^ (note the longer time interval displayed) . 



D. Perfect Homodyne Detection 

In homodyne detection the LO is so strong that indi- 
vidual photons cannot be resolved, and the measurement 
result is a current with white noise £(i). For perfect ho- 
modyne detection the quantum trajectory for the system 
is generated by the SME 



dpi = dt {£ 



£(i)H[e- i$ c 



]}P, 



(2.9) 



which results from the SME in Sec. II C 2 of the preceding 
paper with rj — 1. Here $ is the phase of the LO. 

Homodyne x detection ($ = 0) corresponds to an 'un- 
sharp' measurement of the x quadrature. With a a x 
Hamiltonian, only the measurement causes x to become 
non-zero. Thus, x tends to be projected into an x eigen- 
state, with motion between the two eigenstates on the 
time scale of Y. This behaviour was first noted in Ref. [Q. 
A typical trajectory for SI = 10r is shown in Fig. [|. The 
SI oscillations of y and z are n on-m aximal, and are noisy 
due to the white noise in Eq. (2.9). 

In homodyne y detection (<f> = — 7r/2), the y quadra- 
ture is measured, which pushes the TLA state closer to 
the y eigenstates. Unlike for the x measurement though, 
the SI dynamics quickly rotate y away from the eigen- 
state. as seen in Fig. |J As there is no measurement or 
driving of the x quadrature, it remains strictly zero, so 
the oscillations of y and z are maximal. 
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FIG. 5: A typical trajectory for homodyne y detection with 

n = tor. 



III. LEVELS OF CONDITIONING 

To explicitly show how realistic detection conditions 
the state of the TLA differently from perfect detection, it 
is useful to re-introduce the hypothetical observer of the 
preceding paper (here known as the 'perfect' observer). 
This observer is able to perfectly monitor the system out- 
put as it is absorbed by the realistic detector. A further 
observer is introduced who is able to access only the mi- 
croscopic states of the detector. In the case of the APD 
this corresponds to knowledge of whether the state is 0, 1 
or 2. For the PR, this 'intermediate' observer knows the 
voltage across the capacitor. Hence this observer avoids 
some, but not all, of the detector imperfections. Since 
the detector is being treated classically, the perfect and 
intermediate observers do not affect the average detector 
state because there are no detector superpositions to col- 
lapse. Hence they are consistent with the final 'realistic' 
observer who has access only to the realistic experimen- 
tal records. In the notation of the preceding paper, these 
records are dAf(t) and V(t) for the APD and PR, respec- 
tively. 

Although the measurement records of the three ob- 
servers are obviously different, they can be generated 
from the same experimental run. The observers are mon- 
itoring the TLA at the same time, but have varying lev- 
els of access to the pool of data that could optimally 
be gained from the behaviour of detector. This means 
that correlations between the records exist that will be 
revealed in the evolution of the TLA that they prescribe. 
Through this analysis an understanding of how the spe- 
cific features of the trajectories relate to each detector 
imperfection can be obtained. 

At a particular time each of the three observers will at- 
tach, in general, a different state matrix, p, to the TLA. 
If p represents the state of the TLA in an ontological 
sense then there is clearly a paradox — how could the 



atom have more than one state at the one time? The 
existence of the three observers is, in the context of the 
described experiment, not practical, but it is not difficult 
to think of a realisable experiment in which one observer 
must determine p from a 'filtered' set of measurement 
results, while another has access to the complete set of 
results. Another trivial example is where different ob- 
servers have different knowledge of the initial state of 
the system. However, in this latter case the observers, if 
they had access to the same measurement results, would 
eventually agree on the system state. That is, for a suffi- 
ciently long measurement record, the state is dependent 
upon the record only and not the initial condition. This 
has been shown by Doherty and co-workers || and also 
applies for inefficient detection. 

In Ref. the authors take the view that the state is 
the observer's best description of the system, taking into 
account initial knowledge and the measurement record. 
That is, p represents the state of the system in an episte- 
mological sense. This resolves the paradox raised in the 
previous paragraph, but does not rule out the possibility 
that there also exists a state of the system in an onto- 
logical sense. In particular, the state that the perfect 
observer assigns to the system is unique because it is a 
pure state, and no other pure state can be consistently 
assigned to the system || . A state vector is the best pos- 
sible description of the system (hJ , and for this reason it 
seems that no harm can be done in saying that it repre- 
sents the 'real' state of the quantum system. The states 
assigned by the intermediate and realistic observers can 
only be interpreted as states of knowledge. 

Before leaving this discussion, it is worth saying a few 
words about consistency. It was already noted above 
that although the state matrices assigned by the three 
observers are different, they must be consistent. That is, 
the intermediate observer's state must contain the perfect 
observer's state, and the realistic observer's state must 
contain the intermediate observer's state. When we say 
the state of observer B (pb) contains the state of observer 
A (pa), we can write this as 



which means 



PA C p B , 



3e > : pb - epA > 0, 



(3.1) 



(3.2) 



and could be read as "A knows all that B knows" . For 
example, if the perfect observer knows that the TLA is in 
the ground state, then the realistic observer would have 
to say that the TLA is in a mixture of the ground state 
and some other (not necessarily orthogonal) state. Note 
that this is stronger than the condition for two states to 
be mutually compatible, published recently in Ref. [p|, 
which is simply that they both contain a common pure 
state. 



5 



IV. CONDITIONED DYNAMICS FOR PHOTON 
COUNTING 

In this section we will give quantum trajectories for 
the TLA with output detected using an APD. We con- 
sider both direct detection and adaptive detection using 
a weak LO. We show the system state conditioned on 
three different sets of measurement results, correspond- 
ing to the three observers discussed above. 

The perfect observer's record contains the times at 
which photons arri ve a t the APD. The evolution is there- 
fore given by Eq. ( |2.6| ), with dN(t) known. 

The intermediate observer's record consists of the 
times at which the various APD transitions occur. Thus 
the times at which charged pairs are created, avalanches 
occur and the APD resets are known. When a photon 
from the quantum system is detected this observer's state 
will jump immediately, since the — > 1 transition is mon- 
itored, instead of displaying the delayed jump that is a 
characteristic of realistic detection. 

As a consequence, the intermediate observer has ac- 
cess to results that allow the detector to be described by 
only two states. The evolution of the TLA is the same as 
that which would be the case for a device that has a zero 
response time and a dead time (T dd ) equal to the ran- 
dom response time (r r ) plus the explicit dead time (r dd ). 
The equations that are used to evolve the state matrix 
are similar to the realistic quantum trajectories for the 
APD with response rate 7 r — » oo. These equations com- 
prise the superoperator Kushner-Stratonovich equation 
(SKSE) for the system, which are given in Sec. Ill C of 
the preceding paper, but with T dd replacing the deter- 
ministic r dd . For each avalanche, T dd could be chosen by 
generating a random number R between and 1. Then, 
based on Poissonian statistics, 



r dd = T r + r dd 

= -ln(i?)/ 7r + r dd . 



(4.1) 
(4.2) 



Other changes to the 7 r — * oo equations include re- 
placing dAf(t) by dN cpc (t), which equals one for the time 
intervals in which a charged pair is created (epe), and 
replacing the label of state 2 by dd. This latter change 
indicates that the APD is effectively dead whenever it is 
in this state. The modified equations are given for clarity 



dp 
dpdd 



dt (£ - 7dk - nJ[c + p} + E[dN cpc {t)}) Po 
-dN cpc (t)p + dN cpc {t - r dd )p dd , (4.3) 
dtCp M - dN cpc (t - T dd )p dd 

(■qj[c + fj] +7dk)Po 



+dN cpc (t) 



Tt[(rjJ[c + p] +7 d k)po] 



(4.4) 



with E[dN cpc (t)] = dtTr[(j 1 J[c + p}+ ldk )po}. The super- 
operator C contains the Hamiltonian evolution as well 
as the coupling to the environment, T>[c\. The interme- 
diate observer cannot distinguish between charged pairs 
created by photon absorption or dark counts and is still 
subject to an APD inefficiency. 



The realistic observer has access only to the times at 
which avalanches occur, dftf(t) (and is able to infer the re- 
setting times) . This observer uses the full SKSE to evolve 
his/her state, which we restate here for the reader's con- 
venience 



dpo 
dpi 

dp2 



dt (C - 7dk - riJ[c + p]+ 7 r Tr[pi]) p 
-dJV(t)p + dJV(t-T dd )p2 (4.5) 
dt [(£ - 7r + 7 r Tr[pi]) p x + {r\J\c + p] + 7 dk ) po] 
-dN{t) Pl (4.6) 



dtCp 2 + dM (t) 



Pi 



Tr[pi] 



dM (t - r dd ) p 2 . (4.7) 



As mentioned previously, there will be correlations 
between the three measurement records dN (t) , dN cpc (t) 
and dj\f(t). Numerical simulation will incorporate these 
in a way that will now be discussed. 

At each jump of the perfect trajectory a check is made 
to see whether the detector is in the ready state or not. It 
will be in the ready state if Tr[po] = 1— Tr[p dd ] = 1 for the 
state of the intermediate observer. If it is then a random 
number (/?') between and 1 is chosen to see if the pho- 
ton is actually detected (R' < -q) or not (R 1 > 77). If the 
detection is registered then the stochastic response time, 
r r , is determined accordi ng t o Poissonian statistics, in 



the way indicated in Eqs. ( |4.l| )-(4.2) and the transitions 
that the intermediate observer (dN cpc ) and the realistic 
observer (dM) will register are correlated by 



dN cpc (t) = dM(t + r r ) = dN(t). 



(4.8) 



■ From Eq. (4 



the reader is reminded that an avalanche 
occurs a random time r r after the detected TLA decay. 
The resetting of the APD occurs at the same time for the 
non-perfect trajectories as these observers are monitoring 
(with differing levels of expertise) the same device. If a 
perfect trajectory jump occurs and the APD is not in the 
ready state then no change to the detector takes place. 
A brief discussion of the method of integration of the 
required differential equations for realistic detection will 



be given in Sec. VII 



The x and y quadratures, inversion (z) and purity of 
the TLA are shown in plots of the trajectories. The prob- 
ability of the detector being in each of its states is also 
indicated. An uncertainty in the state of the detector 
exists in the case of realistic detection where the — * 1 
transition is not monitored. 



A. Parameter Values for the APD 

In order that the simulations we perform be relevant 
it is important that we choose a quantum system and a 
detector that are realistic. The quantum system that we 
are considering is a TLA that has an effective damping 
rate given by T — 4g 2 /K, where g is the TLA-cavity cou- 
pling strength and k is the decay rate of the cavity. As 
a guide to the values of these parameters that are real- 
istically obtainable we use the experiment of Turchette, 
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FIG. 6: A typical trajectory for direct detection. The states 
attached to the TLA are plotted for a perfect observer (dotted 
lines), an intermediate observer who can monitor the detector 
transitions (dashed line) and a realistic observer who knows 
only when the avalanche reaches a threshold (solid line). In 
plot (E) the occupation probabilities for the various detector 
states are shown for the above trajectory, with the same leg- 
end applying. The closeness of the lines to the states and 
1 indicates the relative probability of occupation. The real- 
istic parameters used for this photodetection trajectory are 
77 = 80%, and (in terms of the TLA decay rate F = 300Ms _1 ) 
7 r = 7F, r dd = 2r _1 , 7dk = 5 x 10 _6 r and Q = WF. 



Thompson and Kimble p), who have T = 133Ms -1 . We 
choose a value of 300Ms , which is of the same order. 

We must also choose values for 7,, rdd, 7dk and 77. Re- 
alistic values for APDs used in quantum optics laborato- 
ries are Jl|, 7r = 7T, r dd = 2T' 1 , 7dk = 5 x l(T 6 r and 
1] = 0.8. Note that the dark count is negligible. With 
these parameters it will be seen that realistic quantum 
trajectories differ significantly from perfect detection tra- 
jectories. 



B. Direct Detection 

Trajectories for the three observers are given in Fig. ^ 
for direct detection (fj, = 0). The dotted line is for the 
perfect observer trajectory, the dashed line for the inter- 
mediate observer and the solid line for the realistic ob- 
server. The perfect trajectory consists of jumps, which 
mean that the TLA must have decayed to the ground 
state, and no-jump evolution, which causes the TLA 
state to oscillate coherently between the ground and ex- 
cited states. 

When the atom jumps (here for the sake of simplicity, 
we will talk about the perfect trajectory as if that is 
what the TLA is actually doing) then the detector will 
only respond if it is in the ready state and if the photon 
absorption leads to the creation of a charged pair, which 
occurs a fraction 77 of the time. At the time of the first 



TLA jump, the detector is almost certainly in the ready 
state (dark counts are negligible). From y and z in Fig. g 
(B) and (C) we see that the intermediate observer's state 
responds immediately, while the realistic observer has to 
wait until an avalanche builds up before registering the 
detection. By this time, the TLA would no longer be in 
the ground state, so the jump does not take 2 to —1 (or 
y to 0). In fact, the first avalanche actually increases z 
for the realistic observer, contrary to naive expectations. 
This occurs because the response rate (7,. = 7) is smaller 
than the driving strength (SI = 10) for this trajectory. It 
is likely that the TLA will have been driven so that it is 
closer to the excited state than the ground state in the 
time taken for an avalanche to build up. 

After a jump has occurred, the trajectories separate 
significantly. This is because the detector is 'dead' and 
provides no information about the TLA. The evolution 
during this time for the two non-perfect trajectories is 
via the unconditioned ME. The detector can register no 
more photons until it is restored to the ready state, a 
time Tdd after the avalanche has reached the threshold 
level. 

It can be seen from Fig. ^ (D) that when a jump occurs, 
the intermediate observer's state becomes pure (the im- 
purity introduced by the dark counts is negligible). Ob- 
viously, this is because it is known that the TLA is the 
ground state, which is itself a pure state. By contrast, the 
realistic observer's state has a purity of approximately 
0.7 after an avalanche. Apart from the first avalanche, 
this represents a purification. It arises because (incom- 
plete) information has been gained about the time that 
the TLA most recently decayed into the ground state. 

As the response rate 7 r of the detector is large com- 
pared to the photon flux of the TLA (sw T/2 under perfect 
detection) the occupation probability of the ready state 
is close to unity for times after the detector has been re- 
stored to the ready state (see Fig. |6| (E)). That is, the 
transition rate out of state 1 is considerably larger than 
the transition rate into this state. Another feature of 
the trajectories is that the x quadrature remains at zero 
for all time, since the Hamiltonian is proportional to a x . 
Further details and the physical parameters used for the 
detector are included in the figure caption. 

The method of simulating the three trajectories for 
direct detection was to firstly unravel the ME for the 
TLA according to perfect direct detection as in Eq. (tOj). 
Thus, a string of jump times at which dN(t) = 1 was ob- 
tained, along with the trajectory for the perfect observer. 
The method of correlating the other two trajectories has 
already been given in Sec. IV. 



C. Adaptive Detection 

To show how a realistic APD performs when monitor- 
ing the evolution of a TLA in a different way, we cons ider 
the adaptive detection scheme described in Sec. II C . In 



realistic adaptive detection the LO amplitude will be 
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FIG. 7: A typical trajectory for adaptive detection. The sim- 
ulation parameters and legend are as for Fig. pi 



flipped once the avalanche reaches threshold, as opposed 
to the case of perfect detection for which it is flipped 
at an earlier time corresponding to when the TLA de- 
cays. Thus, even the perfect trajectory will not exhibit 
two-state jumping if we flip the LO at the same delayed 
time for all the trajectories. This is done in order to 
stay as close as possible to how the experiment would be 
performed. 

Trajectories for the same three observers as in direct 
detection are shown in Fig. [?]. Looking at the x quadra- 
ture, initially the three states head towards a particular 
eigenstate of the no-jump operator, then a jump occurs 
in the perfect and intermediate trajectories. A short time 
after this the realistic trajectory 'avalanches' and closely 
joins the other two. At this time the LO amplitude is 
flipped. 

The TLA (in perfect detection) then closes in on the 
second no-jump evolution eigenstate, while the detector 
is 'dead'. The non-perfect trajectories are evolving via 
the unconditioned ME, away from the eigenstate. The 
perfect trajectory then jumps away from the x = 1 state, 
implying that a another photon has been detected in the 
combined TLA-LO field. If the LO were being controlled 
by the perfect trajectory, then at this point its amplitude 
would be flipped and the TLA would evolve towards the 
new eigenstate based on the new fi. However, the LO is 
not 'flipped' as this photon goes undetected by the APD. 
The eigenstate of the no-jump evolution, therefore, does 
not change and the TLA in the perfect trajectory must 
evolve back towards it. 

A time Tdd after the avalanche, the detector is ready 
again and the non-perfect states head towards that of the 
perfect trajectory. The next decay, which happens to be 
detected by the APD, takes all the states to the other 
side of the Bloch sphere. 

The y quadrature and z in (B) and (C) display sim- 
ilar features when evolving via the ME or the no-jump 



operator. This is because the ME steady state and the 
operator steady state both have y and z close to zero, 
due to the large f2 being used. As time progresses, the 
three observers' states head towards the steady states, 
thus reducing the y, z oscillation amplitude. 

One can see from (D) that the purity decreases sig- 
nificantly when the APD is dead and unconditioned ME 
evolution is occurring. When it is restored to the ready 
state, the purity increases as no matter where the TLA 
is in the Bloch sphere, it will head towards the no-jump 
evolution eigenstate. The average purity is substantially 
higher than for direct detection. This will be discussed 
later, when we consider system averages. Another point 
is that there is less difference between the purity of the 
trajectories of the intermediate and realistic observers. 
This is because the state of the TLA is not as sensitive 
as direct detection to the time of emission. The dy- 
namics that quickly 'swept' the TLA state away from 
the ground state before the realistic trajectory jumped 
are not present in adaptive detection. 

The simulations of the adaptive trajectories are more 
difficult than those of direct detection, as the perfect tra- 
jectory cannot be run first to determine the decay times. 
This is because the evolution of the TLA depends on 
the LO amplitude, which is dependent upon the realis- 
tic trajectory. However, deciding the jump times of the 
perfect trajectory from those of the realistic trajectory is 
difficult. The best way to do the simulations is to run 
them in 'parallel' so that all three trajectories are being 
evolved at the same time. The perfect trajectory evolves 
under the LO determined by the realistic trajectory, but 
the avalanches of the realistic trajectory are determined 
by the jumps of the perfect trajectory in the same man- 
ner as direct detection. When the APD avalanches the 
LO is inverted for all three trajectories. 



D. Average Conditional Purity 

So far we have shown only typical features of realistic 
detection and how they differ from perfect detection. In 
this section we quantify this difference by considering the 
steady state average purity of the conditional TLA state. 
This is defined as 



p = lim E{Tr[p c 2 (i)]}, 



(4.9) 



where the c subscript is included to emphasize that p c is 
a conditional state. Because purity is a nonlinear func- 



tion of p (2.2), the steady state ensemble average of the 



conditional purity is not the same as the pu rity of the 
steady state ensemble average, which is Eq. (2J5). The 
quality factor p will be less than one, but always greater 
than 1/2. The effect of the detector imperfections can 
be seen by comparing p for the different measurement 
schemes and for a range of driving strengths, f2. 

For small f2, even the unconditioned (without mea- 
surement) stationary purity pme of the TLA approaches 
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FIG. 8: The (scaled and unsealed) ensemble averaged con- 
ditional purity for realistic direct and adaptive detection as 
a function of the TLA driving strength, £1 (in units of F). 
Detector parameters are as for Fig. 0. In this, and subse- 
quent similar plots, the unevenness in the lines is due to, and 
indicates the size of, statistical error. 



unity. To distinguish better detection schemes in this 
limit it is useful to define a scaled purity (between and 
1) that measures how much improvement measurement 
produces: 



Scaled p - 



P- Pme 
1 - Pme 



(4.10) 



For perfect detection the conditional purity is always 1, 
so the scaled p equals 1 also. For no measurement (master 
equation evolution), the scaled p equals 0. 

In Fig. ^| we plot the ensemble averaged purity (and 
scaled purity) as a function of the TLA driving strength, 
fl, for realistic detection, direct and adaptive. The trends 
in this figure result from two main imperfections, the 
dead time and the response time. 

First consider the in put p hoto n flux in each of the two 
schemes, given in Eqs. (2.7)— (2.8) . At low O, direct detec- 



tion produces purer conditioned states. This is because 
in this regime /direct < /adapt., and direct detection is 
therefore less subject to having detections 'missed' due 
to the dead time Tdd- As /direct increases with ft, di- 
rect detection becomes worse than adaptive detection as 

/direct Surpasses /adapt.- 

Note, however, that /direct saturates at T/2, but the 
scaled purity for direct detection continues to decrease 
as ft increases. This is due to the finite detector band- 
width. Under direct detection, the conditioned state of 
the TLA Rabi cycles at rate fi (see Fig. ||). The response 
time 7" 1 can be thought of as the uncertainty in the time 
of photoemission from the TLA, given an avalanche. For 
£1 ~ 7r, the TLA state conditioned on a detection is 
'smeared out' by the consequent uncertainty in how far 
the Rabi cycling has taken the atom from the ground 



state since its emission. By contrast, under adaptive de- 
tection, the scaled purity asymptotes to a nonzero value 
because the conditioned dynamics are governed by T, not 
ft || (see Fig. |?]) . Thus the adaptive purity levels off for 
high Cl. From Fig. ||[ the difference in the purity of the 
trajectories shown for direct and adaptive detection is 
quite large for f2 = 10, as expected from the trajectories 
shown in Figs. || and 0. 



V. CONDITIONED DYNAMICS FOR 
HOMODYNE DETECTION 

For the PR, three different levels of 'realism' will be 
considered, as for the APD. The first observer once again 
has the 'perfect' measurement record. This record is the 
photon flux incident upon the p-i-n photodiode. The 
evolution of this observer's state is via Eq. (2.9). 

The intermediate observer has detailed access to the 
circuit containing the transimpedance amplifier and is 
able to determine / and V (see Fig. 4 in the preceding 
paper). This means, however, that this observer is still 
subject to the diode inefficiency. The current / is given 
by 



I = e ^/p7^ [ V (e-^c + e*V) + y/rj?(t)] , (5.1) 

where the parameters are explained in Sec. IV of the 
preceding paper. The SME for conditioning upon this 
current is 



dpi = dt{C + y /rj^(t)'H[e- i *c]} p 



(5.2) 



The white n oise £'(i) in the above equations is related to 
in Eq. (|H) by 



(5.3) 



This is due to the noise arising from two independent 
sources: the Poisson statistics of the LO and the vacuum 
noise introduced by the inefficiency of the photodiode. 

To relate this to the realistic observer, we use for the 
capacitor voltage 



V' 



V 
RC 



I 

C' 



(5.4) 



or, in terms of the scaled voltage, 



^(e-^c + e^+V^*) 



(5.5) 



Here, V represents the unsealed true capacitor voltage. 
A prime is used to distinguish it from the argument 
(dummy variable) of the probability distribution used by 
the realistic observer. The expectation value in Eq. ( |5.5| ) 
is based on the state matrix obtained through perfect 
detection. 
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The realistic observer will have access to the output 
voltage V measured in the laboratory. The evolution of 
this observer's state is via the SKSE for p(v): 



d Pv (v) = dt[C + ^L^ +1 -v)p(v) 



dv 

+VjdW 3 (t) (v 



(v))p(v). 



(5.6) 



For simulating the realistic trajectory alone, dWj(t) 
would be chosen to be an infinitesimal Wiener increment. 
However, to make it consistent with the perfect and in- 
termediate observers, we must use the following relation 
from the preceding paper. 




V - (v) 



(5.7) 



Ak B T \ 
— dt'y (v) 

^7 dW 3 (t) + dty (v' - (v)) 



dWj(t) 



dt 



(5.8) 
(5.9) 



Here dWj(t) is due to the 'real' Johnson noise, and is 
generated as an infinitesimal Wiener increment. Thus 
v' correlates the equation for the realistic observer with 
those for the other two. 

Simulation of the three correlated trajectories for ho- 
modyne detection was done in parallel, similarly to adap- 
tive detection. This avoids the necessity of storing the 



photon flux and (e **c 



) from the perfect tra- 



jectory at every time step, for later use in the real istic 
trajectories. Once again the reader should see Sec. VII 
for some of the computer programming details. 



A. Parameter Values for the Photoreceiver 

The same quantum system is monitored as in the case 
of photon counting (with T = 300Ms _1 ). For the PR, we 
chose values 7 = 450Ms~\ N = 0.1 and 77 = 0.98. These 
are reasonable values for detectors in quantum optics lab- 
oratories . Note that the efficiencies for the photodi- 
odes of PRs are much higher than those of APDs. This 
is because there are various difficulties in ensuring that 
photon absorptions lead to avalanches. Merely sweeping 
the single charged-pair out of the depletion region is an 
easier task. 

There is generally a trade-off between bandwidth 7 
and the dimensionless noise level N jL2|. We have cho- 
sen a relatively low noise level, and consequently a band- 
width below the maximum available. This noise level can 
be related to the (more usually quoted) noise equivalent 
power (NEP) as follows. Consider for specificity the PR 
model#2007 found in the New Focus catalogue p|. This 
model has a (NEP) of ~ 3pW/VHz (but a bandwidth of 



See attached file HomX.jpg 

FIG. 9: A typical trajectory for homodyne x detection. The 
states attached to the TLA are plotted for a perfect observer 
(dotted lines), an intermediate observer who knows the ca- 
pacitor voltage V (dashed line) and a realistic observer who 
knows only the output voltage V (solid line). Note that at 
some times the dotted and dashed lines are overlapping. In 
plot (E) the occupation probabilities for the scaled capacitor 
voltage (v) are given in a grey scale plot. Darker voltages 
are more likely. The PR parameters are N — 0.1, r\ = 98%, 
7 = 1.5F. System parameters are as for Fig. ^| 



only 790kHz). The NEP can be interpreted as the extra 
optical power that would need to be injected into the re- 
ceiver to simulate noise. To obtain N from the NEP, we 
must compare it to the LO shot noise that is incident in 
homodyne detection. If the LO has transmitted power of 
P then in a time interval St the size of the photon number 
fluctuation will be ~ y/PSt/hujQ . The NEP, on the other 

hand, will produce about (NEPfiwo)v / ^ photons. Thus, 
N = NEP / \/PKJo . For this PR working close to satura- 
tion, P S3 0.5mW 0, so that N S3 0.1 for aA = 780nm 
LO. 



B. Homodyne x Detection 



Trajectories for homodyne detection of the x quadra- 
ture ($ = 0) are given in Fig. |9|. The line-type allocations 
are as for photodetection: dotted is perfect, dashed is the 
intermediate observer and solid is the realistic trajectory. 
The perfect trajectory has been described in Sec. II D. 



Because the photodiodes of PRs have a high efficiency 
(taken to be r\ — 0.98), the trajectories for perfect and 
intermediate detection are very close, especially in y and 
z. The Johnson noise and response time of the circuit 
make it more difficult for the realistic observer to follow 
the TLA state, although a similar trajectory is still ob- 
tained. The smaller detail of the perfect trajectory is 
mostly lost as the realistic observer has trouble identi- 
fying LO fluctuations that have been filtered and then 
obscured by Johnson noise. Another feature is that the 
realistic trajectory never gets as close as the perfect tra- 
jectory to being in an y or z eigenstate, reflecting the 
mixed nature of the state. 

^From the final subplot (E) of Fig. 0, it can be seen that 
to some extent the distribution for the scaled capacitor 
voltage follows the value of x, as one would expect since 
the current is proportional to x. Note that the purity dips 
whenever there is a large fluctuation in the x trajectory. 
This is indicative of quicker evolution being more difficult 
to follow. 
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See attached file HomY.jpg 

FIG. 10: A typical trajectory for homodyne y detection. Leg 
end and parameters are as for Fig. H 



C. Homodyne y Detection 

The trajectories for homodyne detection of the y 
quadrature ($ = — tt/2) are given in Fig. [l0| Once again 
Sec. II D should be referenced for brief comments on the 
perfect trajectory. For homodyne y measurement the 
trajectory associated with inefficient detection is close to 
that of perfect detection, as it was for the homodyne 
x measurement. The realistic trajectory is a reasonable 
approximation to the general shape of the perfect trajec- 
tory, although the y quadrature is not being as closely 
followed as the x quadrature was for the x measurement. 
Note that the amplitude of oscillation of y and z is re- 
duced for realistic detection. 



The distribution for the detector state is influenced by 
y, but oscillations are barely visible. The purity in (D) 
is lower than the purity for homodyne x measurement, 
due the faster evolution of the TLA being more difficult 
to follow. This will now be discussed in more detail. 



D. Average Conditional Purity 

In this section we investigate the long-time ensemble 
averaged conditional purity p as a function of the driving 
strength for the PR. Homodyne measurement of the x 
and y quadratures of the TLA are contrasted. 

The results in Fig. [ll] indicate that as £1 increases, 
homodyne measurement of the y quadrature becomes in- 
creasingly worse at following the evolution of the TLA. 
This is due to the finite bandwidth of the PR in combi- 
nation with the conditional homodyne dynamics in the 
f!>r limit pj . Homodyne y detection produces a con- 
ditional state whose evolution is dominated by fast (f2) 
Rabi cycling (see Fig. |l^). This is because x is strictly 
zero, leaving only the Q dependent y and z in the expres- 
sion for p. 

By contrast, homodyne x detection produces mainly 
slow (r) dynamics, which can still be tracked by the de- 
tector (see Fig. ^). The homodyne x measurement pushes 
x towards the eigenstate (±1), which means that y and z 
must be considerably less than unity as x 1 + y 2 + z 2 < 1. 
The state is, therefore, dependent strongly on x, which is 
devoid of Q oscillation. This explains why increasing f2 
beyond a certain point does not cause further loss of pu- 
rity under x detection. These differences are seen more 
clearly in the scaled purity. 




FIG. 11: The long-time ensemble averaged purity and scaled 
purity is shown for realistic homodyne x and y detection as a 
function of the TLA driving strength, fi. Detector parameters 
are as for Fig. H. 



VI. EFFECTIVE PHOTORECEIVER 
BANDWIDTH 

In Sec. IV B of the preceding paper we presented the 
argument that for small electronic noise N <C 1, the ef- 
fective bandwidth B of a PR is given not by 7 = 1/ RC, 
but rather by 



B 



7 



1-N 



N 



7_ 

N 



(6.1) 



The meaning of B is that we expect the realistic trajecto- 
ries from the PR to be unable to track system dynamics 
which have a rate much larger than B. In this section we 
investigate this claim by studying the realistic quantum 
trajectories for two very different systems. The first is 
the TLA we have used as our model system so far. The 
second is the DPOBT. 



A. Two-Level Atom 

To test the prediction of the effective bandwidth the 
ensemble purity of the TLA was calculated for a range of 
7 and N, while maintaining the proposed quality indica- 
tor of the PR, 7 VI — -W /y/N, as a constant. Our theory 
predicts that the purity will remain constant. The results 
for homodyne x measurement are contained in Fig. [l^ 
and show that the purity is indeed relatively flat, when 
it is considered that N is varying by almost two orders 
of magnitude. 

There is nevertheless a slight downward trend in the 
graph in Fig. [l^. This is most likely due to the presence 
of increased noise as 7 (and N) increases. As explained in 
the preceding paper, our argument only makes sense for 
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FIG. 12: The ensemble averaged purity as a function of 7 for 
realistic homodyne x detection when "fy/1 — N / %/N is kept 
constant at a value of 20. The efficiency is r\ = 98% and 
the driving strength of the TLA is Q = 30. Error bars are 
included to make it clear that the slight downward slope is 
not due to fluctuations. 



small N. In the limit of 7 — > 00 we are returning to the 
situation of adding noise only (discussed in the preceding 
paper) , which is equivalently described by an inefficiency 
of 77/(1 + N). Thus even if B — > 00 as 7 — > 00, a finite 
noise N would still reduce the system purity. 

If B could be kept constant with N decreasing to zero 
then we would expect our prediction to become exact. 
Unfortunately this is a difficult regime to test numerically 
as it is actually the ratio of 7 /AT not j/^/N that appears 
in Eq. ( pOp . With 7, TV — > 0, but B constant, j/N -> 00. 
The time step involved in the simulation decreases and 
the number of iterations of Eq. ( |5.6| ) increases, causing 
the com putational intensity to become prohibitive (see 
Sec. |vij ). 

Further support for identifying B with the effective 
bandwidth is found in Fig. [llj, where purity and scaled 
purity are plotted against £1. It is expected that approxi- 
mately half the total loss of purity will have occurred once 
£1 ~ B. That is, once the frequency of the monitored sig- 
nal becomes equal to the effective inverse response time 
of the receiver. With the parameters of Fig. |ll| we have 
7 VI — N /yN = 4.5, which is in approximate agreement 
with the value of fi such that p = (p max + p m in)/2. Only 
approximate agreement is expected because of the com- 
plicating issue of noise, as discussed previously. 



difficult to do simulations in a regime of theoretical in- 
terest, 7, N — ► such that B is finite. For this reason, we 
now turn to a simpler quantum system, the degenerate 
parametric oscillator below threshold. 

The DPOBT system consists of a damped single-mode 
optical cavity containing a x^ nonlinear crystal which 
is pumped by a (classically described) laser at twice the 
cavity frequency. We take the intensity damping rate to 
be unity and the parametric driving strength, x, to be of 
modulus less than unity. This leads to squeezing in the 
y quadrature, — i(a — a'), if x > and squeezing in the 
x quadrature, a + a', if x < 0. The system obeys the 
following ME 



dp = dt (- 



ixla 2 



p] + V[a]p) 



(6.2) 



when no measurement is performed. 

We will limit ourselves to considering realistic homo- 
dyne measurement of the x quadrature. Using Eq. (5.6) 
with c = a. the SKSE is 



dp v {v) 



dt 



jXH [a 2 



,t21 



■V[a] 




dv V n [ ap(u) + p(u)at ] 

V7dWj(*) (v-{v))p(v). 



(6.3) 



The superoperator can be removed by converting this 
into a Kushner-Stratonovich equation (KSE) for the 
probability distribution, P(x,v), for x and v. This is 
done using the standard procedure for Wigner functions 
presented in Ref. H). The distribution is defined by 



P(x,v) = (x\p(v)\x) 



(6.4) 



with \x) being an eigenstate of the x quadrature oper- 
ator. Due to the linear nature of the optical cavity the 
measurement of the x quadrature does not disturb the 
distribution for the y quadrature after initial transients 
have died away. Thus, all statistics concerning the y 
quadrature can be obtained from the unconditioned ME. 
In particular, the variance of the y quadrature, Ay, in 
the steady state of the ME is 



1 



(1 + x) 2 



(6.5) 



B. Degenerate Parametric Oscillator 

Although the TLA is a small system, amenable to nu- 
merical simulations, it is usually not possible to find ana- 
lytical solutions for its quantum trajectories. This is true 
of perfect detection schemes, let alone realistic detection 
schemes. We have seen above that even numerically it is 



Conversion of Eq. (3.3) gives 
dP v (x, v) 



dt 



k— x - — 

dx 2 dx 2 



7 d 2 




P(x,v) 



^dW 3 (t)(v~(v))P(x,v), (6.6) 
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where, k = ^(1 — x), and the expectation value is found 
from 



= J dxd vvP( X ,v). 



(6.7) 



The damping term has turned into drift and diffusion in 
x, while the parametric driving has become a drift term. 
It should be noted that in this section we are using x 
to represent the possible values that the x quadrature 
can take, rather than being the mean of the quadrature, 
which will be denoted by (x). 

Despite its nonlinear and stochastic nature, a SKE of 



the form of Eq. (6.6) has analytical long-time solutions, a 



fact which is at the heart of modern engineering control 
techniques These solutions are Gaussians, and a 

closed set of equations of motion exist for the conditioned 
mean vector and covariance matrix for x and v. The 
equation for the covariance matrix does not depend on 
the mean vector, and is also deterministic, thus having 
a steady state solution. We denote the three elements of 
the covariance matrix 



A, = (x 2 )-(xf, 

A„ - (v 2 )-(v)\ 
A xv = (xv) - (x) (v) 



(6.8) 

(6.9) 
(6.10) 



It is important to remember the "Ito correction" [p0| in 
calculating the equation of motion for these quantities. 
For example, 



dA x = d (x 2 ) -2{x)d (x) -d(x)d (x) 



(6.11) 



Using integration by parts and the vanishing of P(x, v) 
at infinity yields 




+ + s/ldWj(t)A v 



d(x) = ~k(x)dt + y^dW]{t)A 
d{v) 

dA x 
dA v 



(6.12) 



(-2kA x + 1 - 1 A 2 xv )dt 

7 o nv 



(6.13) 
(6.14) 



N 




2 1 A V - 1 A 2 V dt (6.15) 



(k + j)A xv + 




l) + jA v A xv 



dt, (6.16) 



which are examples of Kalman filter equations [Q. In 
the long time limit the covariances are constant, and the 
whole distribution is just being shifted with the motion 
of the mean. 

Since the purity of the cavity mode is defined in terms 
of the variances of the x and y quadratures, after initial 
transients have died the purity will be constant in time. 
For Gaussian states with no correlation between the x 



and y quadrature, the purity Tr[p 2 ] is given by |Ti 



1 



P = 



v / A x A y 



(6.17) 



There is an analytical steady state solution of 

Eqs. (6.14)-( |3.16j ), but it is very complex and will not 
be given here. We are really only interested in the purity 
in the limit of small N with B — j/\/~N fixed. If B is in- 
deed the effective bandwidth, as we have argued, then in 
the limit N — ► with B fixed, the purity should depend 
on B only, not 7 or JV. 



To show this we will examine Eqs. (6.14)— ( 3.16) more 
closely. The first step is to determ ine h o w , A„ , A^,, 
scale as 7, N — > 0. Re- writing E qs. ( | 6. 14 )— ( 5. 16 ) in terms 
of B and N, we have from Eq. ( 6.14| ) at steady state 



A., 



1 - bVnaI, 

2k 



(6.18) 



^From the unconditioned ME, A^ would be l/2fc, while 
we expect (and have verified numerically) that a finite B 
will give a finite variation of the purity away from pme- 
For this to be the case it must be true (when 7, N — ► 0) 
that A x < I /2k with 



BVNA 



NAt 



1 



A xv ~ N- 1 / 4 . (6.19) 



Ignoring the 7 term in Eq. (6.16) that is small compared 
to k gives 



1 



BVN A xv NV*y B 



1) 



where we have used Eq. ( 6.1 9| ) . 

New variances of order unity are now defined: 

A v = N 1 / 2 A V 
A xv = N^Axv. 



(6.20) 



(6.21) 
(6.22) 



The variance in x is kept the same , as it is of the order 
unity. Using these in Eqs. ( |6.14 )-( |6~T6 ) and ignoring 
small terms allows these equations at steady state to be 
written as 

-2kA x + l- BAl v = (6.23) 
B - 2yJ~BqA xv — BA 2 V = (6.24) 
-kA xv - \f~Bri(A x - 1) - BA XV A v = 0. (6.25) 

Importantly, these equations, which apply in the 7, N — > 
limit, are only dependent upon the detector parameters 
B and ij. This means that the 7, N dependence of the 
purity can be summarized by the one parameter B = 
7/viV. This confirms the correctness of our argument. 

As Eqs . (|6.23| ) ( 3.25 ) are more simple than 
Eqs. (6.14)-(p.l6[) it would be useful to solve for 



A X: and hence purity. This can be done, using the facts 
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that the variances are positive and that A^Aj, > 1 to 
discard non-physical solutions. The result for purity is 



where 



2fcr?(l - k)B 2 



1/2 



Sy^qB* {BR + k 2 ) - fc 4 - B 2 k 2 (l - f + 



with 



R = ^k 2 + n{l-2k) 

S = k^(2BR+k 2 + B 2 )/T 1 B 3 



(6.26) 



(6.27) 



We can now consider the two obvious limits of B — > 
0,oo (with rj = 1 for simplicity). For B — > the purity 
is, to lowest order in B, 



P = Pme + 



B 2 (l-2fc) 2 (l + fc) 2 
~4 F/ 2 (l - fc) 3 / 2 



As _B — > oo the purity to first order in 1 /B is 



p=l 



1 (l-2fc) s 
LB (1-fc) 



(6.28) 



(6.29) 



That is, for small J5, the purity increases from the un- 
conditioned purity quadratically with the effective band- 
width, while for large B it decreases from unity linearly 
in B~ x . 

As a final point, it is worth noting that the down- 
ward trend observed in the purity in Fig. |l2| for the 
TLA is also bo rne out in the DPOBT. That is, for fixed 
B = 7vl — N /N, the purity still has a slight downward 
trend with N. This was found through numeric investi- 
gation of the complete solution for the purity (which was 
not stated due its complexity). However, the numeric in- 
vestigation also revealed that as N — > the purity curve 
maintains the same slope, going to the limits found ana- 
lytically above. 



VII. NUMERICAL SIMULATIONS TECHNIQUE 

Before concluding, we here comment explicitly upon 
the method of numerical simulation of the SKSE's con- 
tained in this paper. Obtaining realistic trajectories from 
Eqs. (]4 . 5|) — ( 4 . 7 ) and Eq. (5.6) is obviously a non-trivial 
numerical exercise, even when they do not have to be 
correlated with the unrealistic trajectories of the perfect 
and intermediate observers. To set up the problem, the 
Quantum Optics toolbox for Matlab p8[ was used. This 
allowed easy formulation of the required quantum (su- 
per)operator expressions. 

We represented the supersystem state with a long col- 
umn vector and a square matrix, L, was used to evolve 
it. That is, 



{Pi\ 



V J 



= dtL 



V- J 



(7.1) 



Ps = 




(7.2) 



The subscripts e, g represent the ground and excited state 
of the TLA. The integer subscripts s label detector states 
s 6 S. For the APD there are only three detector states, 
while for the PR the scaled capacitor voltage is dis- 
cretized on a grid. The m atri x L ge nerates al l the evolu- 
tion on the RHS of Eqs. (U)-(|j) and Eq. (fl|. 

Because of the stochastic nature of the realistic tra- 
jectories, L cannot in general be formed in its entirety 
at the start of the simulation. In the case of the APD 
one can avoid t his d if ficul ty by using the unnormalized 
versions of Eqs. (U)-(fO). These 

were given in |2J . The 
stochasticity then enters via the comparison of the norm 
of the supersystem state to a random number in order 
to choose the avalanche times. For the PR, the nonlin- 
ear term (the expectation value, (v)) is included in the 
evolution. This, in addition to the noise, c?Wj(t), means 
that the best that can be don e is to create the L for all 
but the last term of Eq. ( |5.6|) . The last term is created 
every time step after the calculation of the expectation 
value. 

The structuring of the numeric solution in this way 
(creating as much as possible of L before the iteration 
process begins) makes it very flexible. For example, a 
change in the nature of the Hamiltonian for the TLA 
could be easily achieved by creating the new superoper- 
ator in the Quantum Optics toolbox and thus obtaining 
a new L. This is in contrast to the technique of deriving 
x s ,y Sl z s equations and then evolving these numerically, 
as a new Hamiltonian would require a new derivation. 

The obvious disadvantage with L is that it will be very 
sparse in general. This is overcome by reducing it to only 
the non-zero elements, with the aid of the find command 
in Matlab. Once the reduced matrix has been found it 
is written to a data file which is used as the input for a 
C++ program. 

The C++ program then integrates the matrix equation 
7.1 by looping through all the non-zero L elements and 
making the appropriate increments. The simple Euler 
technique of integration is u sed, by which is meant that 
the infinitesimal d in Eq. (7.1) is replaced by the small but 
finite 8. This method has well known instabilities, but is 
used by other workers in the field of quantum optics. It 
has the property of being accurate when the solutions do 
not explode, as opposed to more stable methods which 
have less spectacular failures that are more difficult to 
detect 

Some specific details of the PR simulations are as fol- 
lows. A stationary grid of 100 points was used for the 
capacitor voltage distribution. These points were spread 
7 standard deviations of the initial distribution either side 
of the initial mean. The initial distribution was found by 
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solving the Ornstein-Uhlenbeck equation MM 



VIII. CONCLUSION 



(7.3) 



This is derived from the SKSE (5.6) by removing the 



TLA and averaging over the realistic measurement. Ig- 
noring the TLA is reasonable as the field from the TLA 
is of the same order as the vacuum field. This leads to 
a Gaussian solution with a variance of A v u — 1/2N, 
giving a standard deviation of 2.24 for N = 0.1. The u 
subscript is to indicate that this is a variance uncondi- 
tioned on measurement. From Figs. (P)~(|lOl), it can be 
seen that the considered range of voltages was sufficient. 

A moving grid for the voltage distribution was consid- 
ered but not used in the end. This decision can be justi- 
fied by calculating the conditioned variance of the voltage 
distribution in the case of a vacuum input. Unlike the 
calculation of the variance from the Ornstein-Uhlenbeck 
equation, the stochastic measurement term 



VldWj(t) (v-(v))P(v) 



(7.4) 



is now included. The variance goes to a steady state 
value, despite the stochasticity. The derivation of the 
conditioned variance, A^.c, is performed in the same 
manner as in Sec. VI B, and the result is 



A, 



1 

N 



(7.5) 



This is always smaller than the unconditioned variance, 
1/2N. In fact, in the limit N -> 



1 



2VN 



(7.6) 



In this limit the conditioned variance becomes much 
smaller than the unconditioned variance and a moving 
grid would save much numerical computation. This is 
because the number of grid points necessary to describe 
the non-zero probabilities at a particular time is much 
less than those required to describe the movement of the 
distribution over all time. However, for TV = 0.1 we have 
A„ iU = 2.16A t , !C and the saving is not very large. It is 
worth noting that a moving grid would not solve all the 
computational intensity of j/N — > as the problem of 
the decreasing required time step still exists. 

For the PR, time steps of St = 1 x 10 _5 r _1 generally 
proved satisfactory as did ensemble sizes of about 1000 
for forming averages. The time step for the APD can 
be increased by about an order of magnitude as there is 
no white noise. The ensembles actually took the form 
of a collection of samples of the supersystem state in 
one trajectory, taken at times separated by This is 
large compared to the system correlation time |24j . The 
equivalence of this time averaging to a many trajectories 
average has been established by Cresser [pTl . 



In this paper we have applied the theory of realis- 
tic quantum trajectories for photodetection derived in 
the preceding paper |jj to realistic quantum optical sit- 
uations. We investigated two systems. The first was 
a driven two-level system, which could be realized as 
a strongly coupled atom-cavity system, heavily damped 
though one cavity mirror. For this system we solved for 
the trajectories numerically. We looked at features of 
typical stochastic trajectories, as well as a figure-of-merit, 
the ensemble averaged conditional purity, for how well 
the system state is known. We considered four different 
detection schemes: direct and adaptive detection using 
an avalanche photodiodc (APD), and x and y quadra- 
ture homodyne detection using a photoreceiver (PR). 
The second system was a below threshold degenerate op- 
tical parametric oscillator. The conditioned evolution for 
this could be solved analytically for homodyne detection 
using a photoreceiver. 

The first significant result we found is that for realistic 
detector parameters and realistic systems there is a large 
difference between standard (idealized) quantum trajec- 
tories and realistic quantum trajectories. Even with in- 
efficiencies and dead time included in the standard tra- 
jectories, the differences are marked. This is due to the 
effect of finite detector response times. For example, for 
the strongly driven two-level system, the average condi- 
tional purity under realistic direct detection was scarcely 
better than with no detection at all. 

The second significant result we found was the amount 
of difference the detection scheme makes, even using the 
same detector. With the APD, the adaptive scheme gave 
(for most parameter regimes) far higher purity than di- 
rect detection. For the PR, homodyne x detection (with 
a local oscillator in quadrature with the system driving 
field) gave a similarly better result than homodyne y de- 
tection. With a perfect detector the choice of detection 
scheme of course makes no difference as the conditional 
purity would always be unity. 

A final significant result we found was that for homo- 
dyne detection the PR bandwidth 7 is not the relevant 
parameter for determining the rates of system evolution 
that can be tracked. Instead, the effective bandwidth is 
B = j/yfN. Here N is the ratio of electronic noise to 
vacuum noise, which has been assumed small (as required 
in quantum optics experiments). We verified this result 
numerically for the two-level system and analytically for 
the parametric system. We would expect this result to 
be generalizable to other sorts of detection scheme, as 
long as they involve filtering and the addition of white 
noise. Mesoscopic electronics is an obvious example, as 
discussed in the preceding paper. 

All of these results are relevant to the area of quantum 
control. As explained in the introduction, since a condi- 
tional state is a representation of the observer's knowl- 
edge about a system, it is by definition the optimal ba- 
sis for controlling that system. For this to work, the 
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observer must have an accurate model for the relation 
between the available information (the detector output) 
and the quantum system. That is precisely what a realis- 
tic quantum trajectory is. The large differences between 
standard and realistic quantum trajectories noted above 
means that the former would be a poor approximation 
to the latter in a control system. That is, a realistic 
quantum trajectory theory is necessary. The differential 
performance of different measurement schemes may also 



be significant, if one can choose the measurement scheme 
to be used in the control loop. 

Finally, another area where realistic quantum trajecto- 
ries would probably be needed is in the estimation of dy- 
namical parameters for open quantum system from mon- 
itoring the outputs of such systems. This has been inves- 
tigated for perfect detection in Refs. |22|, The loss 
of information resulting from realistic detector imperfec- 
tions will be the subject of future work. 
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